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Abstract 

Background: The neutrinoproduction of photons and pions from nucleons and nuclei is relevant 
to the background analysis in neutrino-oscillation experiments [for example, MiniBooNE; 
A. A. Aquilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev. Lett. 100, 032301 
(2008)]. The production from nucleons and incoherent production with E u ^ 0.5 GeV have 
been studied in [B. D. Serot and X. Zhang, Phys. Rev. C 86, (2012) 015501; and X. Zhang 
and B. D. Serot, Phys. Rev. C 86, (2012) 035502]. 

Purpose: Study coherent productions with E v ^ 0.5 GeV. Also address the contributions of two 
contact terms in Neutral Current (NC) photon production that are partially related to the 
proposed anomalous uj(p), Z boson and photon interactions. 

Methods: We work in the framework of a Lorentz-covariant effective field theory (EFT), which 
contains nucleons, pions, the Delta (1232) (As), isoscalar scalar (a) and vector (uj) fields, 
and isovector vector (p) fields, and incorporates a nonlinear realization of (approximate) 
SU(2) L (g) SU(2) R chiral symmetry. A revised version of the so-called "optimal approxima- 
tion" is applied, where one-nucleon interaction amplitude is factorized out and the medium- 
modifications and pion wave function distortion are included. The calculation is tested 
against the coherent pion photoproduction data. 

Results: The computation shows an agreement with the pion photoproduction data, although pre- 
cisely determining the A modification is entangled with one mentioned contact term. The 
uncertainty in the A modification leads to uncertainties in both pion and photon neutrino- 
productions. In addition, the contact term plays a significant role in NC photon production. 

Conclusions: First, the contact term increases NC photon production by ~ 10% assuming a 
reasonable range of the contact coupling, which however seems not significant enough to 
explain the MiniBooNE excess. A high energy computation is needed to gain a firm conclu- 
sion and will be presented elsewhere. Second, the behavior of coherent neutrinoproductions 
computed here is significantly different from the expectation at high energy by ignoring the 
vector current contribution. 
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I. INTRODUCTION 



This is the third of a series of studies about neutrinoproduction of photons and pions with 
neutrino energy E v ^ 0.5 GeV where the A excitation is important |l|-|3|]. The focus of this 
paper is the coherent production. As we know, in the neutrino-oscillation experiments, for 
example MiniBooNE |4H6|, the photon and pion neutrinoproduction from nuclei and nucleons 
are potential backgrounds. It is still a question whether NC photon production might explain 
the excess events seen at low reconstructed neutrino energies in the MiniBooNE experiment, 
which the MicroBooNE experiment plans to answer 0. In Refs. HHnl, the authors argued 
that the anomalous interaction terms involving u(p), Z boson, and the photon may increase 
neutral current (NC) photon production. So the cross-section calculation for these processes 
becomes necessary. 

In Ref. [H, we introduce the A resonance as a manifest degrees of freedom in a Lorentz- 
covariant effective field theory (EFT) with a nonlinear realization of the SU(2) L <g> SU(2) R 
chiral symmetry 1 . In Ref. [2J, we study both neutrinoproductions from free nucleons and 
calibrate our theory. Because of various symmetries that are built in, the conservation of 
vector current (CVC), and the partial conservation of axial current (PCAC) are satisfied 
automatically, which is crucial for photon production calculation. In Ref. [jjj], we work on the 
incoherent productions from the nucleus. The previous studies show that the contributions 
due to the two terms mentioned above [Io| are tiny in the NC photon production from both 
free nucleons and nuclei. This paper is devoted to the study of coherent productions from 
12 C, which is the major target nucleus in the MiniBooNE's detector, and also to addressing 
the significance of the two mentioned terms. 

Here we apply the so-called "optimal" approximation, in which one-nucleon interaction 
amplitude can be factorized out from the full nuclear matrix element leading to great sim- 
plification of the calculation. Meanwhile both the CVC and the PCAC are preserved. The 
nuclear ground state is calculated by using the mean-field approximation (see Ref. [3] for 
the details). The o ptim al approximation was first illustrated generally for projectile- nucleus 



scattering in Refs. [20|, [21 



It has been applied quite successfully to nucleon-nucleus scatter- 
ing in a relativistic framework 22|-25j . Moreover, a similar approximation has been applied 
in pion -nucleus elastic scattering [26M28 



ion-nucleus elastic scattering 
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and coherent pion photo- and electro-production 
2 . It was realized that the medium- modification of the one-nucleon interaction am- 
plitude plays a key role in the coherent production, which is also included in our revised 
approximation. In the QHD EFT model, baryons interact with each other via. exchanging 
the mesons (in the space-like region). Because these bosonic fields develop finite expectation 
values in the medium, the real part of the baryon self-energy is modified on the mean-field 
level Meanwhile the change of the A width has been studied in the nonrelativistic 

framework both phenomenologicall y |36l . 37| and theoretically (38|, but it is not completed 
in the relativistic framework [3|, |39|-|42| . Here we continue our simple treatment proposed in 
Ref. [3|. Another important factor related to the A is the distorted pion wave function in 
the pion production, which is included here by using the Eikonal approximation. Such an 
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19( , and is often called quanti 



1 The EFT was originally motivated by the nuclear many-body problem 
hadrodynamics or QHD. 

2 Ref. [351 ] pointed out that using different one-nucleon interaction amplitudes that are equivalent on shell 
can lead to quite different results. Here such ambiguity does not exist because we have a unique free 
interaction amplitude. This will be addressed later. 
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effect can be ignored for photon production. Comparing the two may be used to disentangle 
the medium-modification and the pion wave function distortion. 

To benchmark the approximation scheme, we calculate various differential cross sections 
for pion photoproduction. We are able to get an agreement with existing data (431-45 



The approximation is then applied to study the photon and pion neutrinoproduction. Un- 
fortunately, existing neutrino experiments, for example Refs. 0, 47||, do not put a strong 
constraint on pion productions with E u ^ 0.5 GeV, since most of them have only spectrum- 
averaged measurements, and the mean neutrino energy is around 1 GeV. On the theoretical 
side, there are other microscopic calculations on pion productions (48T-58]. In most of them, 



the optimal approximation is in one way or another applied. The A dynamics is taken into 
account by using the nonrelativistic models. The final pion wave function is calculated either 
in the Eikonal approximation or by solving the Schroedinger equation with the pion optical 
potential. The key difference between our work and others is that we work in a Lorentz- 
covariant EFT, which has been applied successfully to nuclear many-body problems and also 
has been calibrated for neutrinoproductions from free nucleons. The medium-modification 
of baryons can be calculated on the mean-field level. We can address the power counting 
of different diagrams in this EFT, although the theory can only be used at the low energy 
region (E v ^ 0.5GeV). More importantly, coherent NC photon production has rarely been 
discussed in the microscopic approach. In addition, there exists a macroscopic approach, 
which treats the nucleus as a whole and makes use of the forward scattering behavior of 
coherent pion production in the high energy scattering. In the forward scattering kinemat- 
ics, PCAC leads to a relation between the pion neutrinoproduction and pion-nucleus elastic 



ie_p 

scattering. This is initiated in Refs. [59|, |60|, used in the NUANCE event generator [6JJ, |62 



and revised recently in Refs. [63H65 



This approach has also been applied to compute coher- 
ent photon neutrinoproduction at the 2 GeV region and beyond [66 1. The NUANCE output 
on coherent pion production shown throughout this paper, with which we will compare our 

results, is obtained from the NUANCE v3 event generator with the calibration applied from 
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the experimental data [47J, [61|, 167 



The paper is organized as follows. In Sec. [TTJ the approximation scheme is discussed, 
and the difference between our approach and others is emphasized. Sec. II III presents our 
results. At first, differential cross sections for pion photoproduction are compared to the 
data. The effect of the two contact terms are discussed. Then, we show the results for pion 
neutrinoproductions. Since there is the uncertainty in our model, results of using different 
parameters are compared. Finally, we focus on NC photon production and discuss the 
relevance of our results to the MiniBooNE low reconstructed energy excess events. The two 
contact terms are again discussed in this context. 



3 The experimental analysis indicates that to be consistent with the coherent NC pion production data, a 
35% reduction needs to be applied to the orginal NUANCE output 



47 



67|. 
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II. APPROXIMATION SCHEME 



A. Kinematics 

The formulas needed for computations are shown here. For the 7r° photoproduction, 



a 



*-! m \n, . (!) 



^7r) 2 E 1 m A 



Here g M and k w are the momentum of the incoming photon and outgoing pion, and q° = J5 7 
is the photon energy in the laboratory frame. Because the nucleus A remains in the ground 
state and is heavy enough to ignore its recoil, we can have = q . The 1/tua term is used 
to properly normalize the quantum state. The definition of transition probability is 



= ^E e 1 c ^9)<A7r(4)|^> 2 > (2) 
m A m A z x _ 

where J% ad is the Electromagnetic current involved in this process, and A, is the photon 
polarization. 

For the neutrinoproduction, 



We define pu and pif as the momenta of incoming and outgoing leptons. q = pu — pif, 
p\ = E u (the incoming lepton energy in the laboratory frame). Here we also have k® = q°. 
The nuclear matrix element is 



|M| 



2 



m A m A 



\Y.^ V u<iGf) 2 l»{p lu p l} )(AMk)\Jl d \A) 2 • (4) 



Here Gp is the Fermi constant; V u d is the u and d quark mixing in the charged current 
(CC), and is 1 in the NC; and l^,(pn,Pif) is the corresponding lepton current. Moreover, 
J£ ad = r EM = Ijg + for the photoproduction, J^ ad = \{V^ + A*), i = ±1 for the CC, 

and J Ld = \(Y^ + A °^) ~ sin 2 Ov>Jem for the NC is the weak mixin g angle). 

For NC photon production, the zero mass of photon should be taken into account in 
Eq. ©, and (A, tt\ J£ ad \A) in Eq. (0} needs to be changed to (A,~/\Jg ad \A). 

B. The optimal approximation 

The current matrix element can be written as 

■ j A dre^-^iJlLM #0) PW, 

-> ■ poo Yl(pJ) ii 

(5) 

^dfe^-^e tJz 2 i^i (J£ a M^,r)) DW. 
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FIG. 1: The Feynman diagrams for pion production. Here C stands for various types of currents 
including vector, axial-vector, and baryon currents. Some diagrams may be zero for some specific 
type of current. See Ref. [2j for the details. 



Only the one-body current contributions are included coherently. We apply the optimal 



approximation to simplify the calculation [29M35 



Z , Pni Z Z 

+ Pp^\Y^{P^ S ^^-^VKaA^)\p,S z ,^—^) . (6) 
Z . Pni Z Z 



Refs. [20J, [21| argued that in the center mass frame of the projectile and the nucleus, the 
nuclear matrix element can be expressed as the product of a proper density and the free 
one-nucleon interaction amplitude calculated in the Breit frame of the projectile and the 
nucleon. Ignoring the recoil of the nucleus leads to Eq. (EJ. For NC photon production, we 
can use Eqs. flSJ) and (EJ) with a proper current inserted. 

The calculation of the one-body current matrix element for both pion and photon pro- 
duction in Eq. has been discussed in Refs. j^, 0] • There are two basic types of Feynman 
diagrams contributing here, as shown in Fig. [1] diagrams with the A [(a) and (b)] and all 
the rest called nonresonant diagrams here. The diagrams for the photon production can be 
viewed as those in Fig. [1] with the final pion line changed to the photon line. The medium 
modification on the one-nucleon interaction amplitude, as introduced in Ref. is based 
on the mean-field approximation. The effective mass is introduced for the baryon to include 
the modification on the real part of its self-energy: 

M* = M - g s ((j)) , (7) 
m* = m — h s ((fi) , (8) 



P°n = Pn + 9v(V°) = VM*z+i? n + g v (V°) , (9) 



p° A = pf + h v (V°) = ^m* 2 +f A + h v (V°) . (10) 

Here g SyV (h s>v ) are the couplings between the scalar and vector mesons and the nucleon (the 
A). Figs. [2] and |3] show the calculated g s {4>) an d g v (V°) m 12 C (we approximate it as a 
spherical nucleus). "Gl" and "G2" label two parameter sets about g s , g v , and others 
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In this paper, we use "Gl" as in Ref. jjjj]. For the A width, we follow Refs. j3, 36, 37]. Above 
the pion threshold, 

Ta = + r sp j 

p{r) 



sp 



V x 



P(0) 



where is the A pion-decay width in the nucleus, and can be found in Refs. [38|, |4l|, |68| . r sp 
is the width of other channels, which has been fitted in Refs. 37|, |56j. We set Vo ~ 80 MeV 
(!, 37, 56]. p(r) is the baryon density at radius r. Below the pion threshold in the photon 
production, 

r A ~ r sp « v x 4^ • 

In the cross channel of the A diagram, we set the width to be zero. 
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FIG. 2: (Color online). The proton and neutron densities in 12 C as calculated in the mean- field 
approximation by using Gl and G2 parameter sets 14]. 



In addition, since pions interact strongly with the nucleus, it is necessary to treat the 
final pion wave function in a realistic way. As shown in Eq. fl5]), the Eikonal approximation 
is used to calculate the distorted wave function |48j, which is labeled as DW, while the PW 
calculation is without such distortion. For NC photon production, we only apply the PW 
calculation. In Eq. (jSJ), U(p(f),z) is the pion polarization insertion in the nuclear medium 
with baryon density p(r), as calculated in the local Fermi gas approximation. Following 



Refs. 



49, 51 



n 



v 



we use the following formula for II in symmetric nuclear matter: 
V 



-4tt kl 

s *1 + Airg'V 



m 



Re^Ao + — HmT^A, 



-v/s-m + 2M-ReS 



AO J 



(11) 
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where g' = 0.63, £a is the A self-energy insertion, and T n is the A pion decay width as 



discussed before. We take the results from Ref. [49] for the Ea and ReEAo (See Refs. [68l. 169 
for the details) 4 . 



C. The approximation used in Ref. 35] 



It is interesting to compare our calculation with that in Ref. [35] where the relativistic 
mean- field theory is also used. Instead of using Eq. (JHJ), the authors there project the one- 
nucleon interaction amplitude to an independent basis, and then convolute the amplitudes 
of each basis with the corresponding current densities calculated in the relativistic mean- 
field theory. Take the proton contributions, for instance. First decompose the free proton 
interaction matrix element: 

(p. ^^\JLM fcOb. = % + K a i a + + ...) Ul , (12) 

and then multiply the amplitude, for example Fy a , with the proton vector current density 
{A\il}p(x)'y a ip p (x)\A) . The sum of different terms' contributions in Eq. ( fl2l) is the proton 
contribution to the nuclear matrix element. For the closed shell nucleus, the only relevant 
amplitudes are F s , F v and F T , because the densities associated with other amplitudes in 
Eq. ( I12p are zero for a spherical nucleus. We have compared the calculations for the pion 
production by using Eq. ( TT2|) with those by using Eq. (JBJ). Fig. H] shows the comparison 
of the total cross section of coherent n° photoproduction: the "current decomposed" uses 



4 We essentially treat Eq. (fTTj) as an analytical expression for the pion optical potential, but it is not used 
to deal with the A modification in the one-nucleon interaction amplitude. 
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FIG. 4: (Color online). The photon energy dependence of the total cross section for coherent 7r° 
photoproduction from 12 C. In both calculations, we set (r s , r v ) = (1, 1), which controls the A 
medium modification. The explanation for different calculations can be found in the text. 

Eq. (TT2T) . while the other uses Eq. (J6j). In the two calculations, we include the same medium- 
modification to the one nucleon matrix element (Ref. [35J uses the free amplitude) and 
the same Eikonal approximation to calculate final pion wave function. For the A medium 
modification, we set (r s , r v ) = (1, 1) [r s = h s /g s , r v = h v /g v ; see Eqs. ((Tj) to ({TO]) ]. Only 
the A diagrams are considered, because including the others needs the extra care of the 
Electromagnetic current conservation in the "current decomposed" calculation. We see the 
difference between the two is small. As we have checked, this is also true for the differential 
cross section and for the cross section of the neutrinoproduction. 



In addition, it was pointed out in Ref. [35( that there exist other amplitudes that have 
the same on-shell behavior in nucleon scattering but give quite different results for nucleus 
scattering through using Eq. (j!2p . There is no such ambiguity in our approach, because we 
have a unique interaction amplitude derived from the QHD EFT Lagrangian. This shows the 
importance of having a consistent framework describing both nucleon and nucleus scattering. 

III. RESULTS 

A. Coherent it photoproduction 

Fig. [5] shows five different calculations for the photon energy dependence of da/dQ n of 
pion photoproduction from 12 C with the pion angle fixed at 6 n = 60° ± 10° (relative to 
the incoming photon direction). All the variables are measured in the laboratory frame of 
the nucleus. The data are from Ref. [431 ] . These calculations include diagrams up to the 
v = 3 order [y = 2 terms do not contribute in this production) [2[. Again in the labellings 
of different curves, (r s , r v ) are defined as (h s /g s , h v /g v ). Since the two couplings are not 
precisely known, we simply show results with three different choices: (0, 0), (1, 1) and 
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FIG. 5: (Color online). The photon energy Ey dependence of da/dQ w for coherent 7r° photopro- 



duction from 12 C. The final pion angle is fixed at 6 n = 60° ± 10° 
calculations can be found in text. The data are from [43(. 



The explanation for different 



(1, 0.9) 5 . As discussed in Refs. 
photon, nucleon, and Z boson ( 
of the photon: 



Or 7T 



, there are two low-energy contact terms involving the 
that contribute at the v = 3 order in NC production 



Cl 

M 2 



N^NTt 



ei 
M 2 



N^a u Nf 



SflU 



Here F^J and f v are related to the photon field, and a u is related to both Z boson and pion 
fields. Interestingly in Ref. [33J, it is shown that c\ term plays a significant role in coherent 
pion photoproduction. Refs. [8Hl0j point out the anomalous interactions of to and p° mesons 
can induce such contact terms at low energy with c\ = 1.5 and e% = 0.8 6 . However, as 
argued in Ref. 0, c\ can also be induced by the off-shell interactions involving the A, which 
leaves its value unfixed. In our calculations shown in Fig. [51 we use c\ = 1.5, except for 
those labeled with il (c% = 3)" where we double C\. (Since the e± term's contribution vanishes 
for an isospin target, we focus on C\ in the following.) 

We can see that the first two calculations with "(0 , 0)" fail to give the right predictions 
around the peak. "(1, l),u = 3 (ci = 3)" and "(1, 0.9), v = 3" give the best predictions. 



In Ref. [33|, it is also noticed that fixing the real part of the A self-energy is correlated with 
c\. However when the photon energy is above 0.3 GeV, all the calculations underestimate 
the cross section. The shapes of different curves are controlled by the nuclear form factor, 
e.g. the Fourier-transformation of nuclear densities [see Eqs. (jSJ) and (jSJ)], and hence they 
are similar. So we expect the underestimation to be generic for all the one-body-current 
calculations. To resolve this issue, two-body currents may need to be considered. In addition, 
around the peak, the (0, 0) result is smaller than the (1, 1) and the (1, 1) smaller than 
the (1, 0.9) [ca = 1.5 in all the three). The same pattern has been found in the incoherent 
productions [3j. It was argued that among the three, the (0, 0) requires the most energy to 



5 It has been shown that (r s , r v ) = (0,0) can not explain the A spin-orbit coupling. Detailed discussions 
about the_two can be found in 



6 In Ref. 



33| . the authors use the w's pion-decay vertex to generate the c\ coupling {c\ is also around 1.5). 
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excite the A, while the (1, 0.9) requires the least. In the coherent production, the nuclear 
form factor makes them even more sensitive to r s and r v . This will also be seen for the total 
cross section of pion neutrinoproduction. 

In Fig. 0, we show our calculations for the scattering angle 9 n dependence of da/dfl n with 
the photon energy fixed at = 0.173, 0.235, and 0.29 GeV. All the variables are measured 



in the laboratory frame of the nucleus. The data are from Refs. [44], |45j. Each plot shows the 



same five calculations as those in Fig. [5J Systematically with c\ = 1.5, the (0, 0) prediction 
is smaller than the (1, 1) and the (1, 1) smaller than the (1, 0.9). In the forward kinematic 
region, i.e. small 9 n , both "(1, 1), v — 3 (ci = 3)" and "(1, 0.9), v = 3" agrees with data for 
the three cases. However for larger 6 n , the calculations fail: for E 1 = 0.235 GeV, the two 
overestimate the cross section when 20° ^ 9^ ^ 60° and underestimate it when 8 n ^ 60°; 
for E 1 = 0.29 GeV, the two give too big results compared to the data when 6 n ^ 40°. 
Nevertheless, we expect our calculations to work better at the higher energy region, because 
the cross section is more dominated by the forward production. 



B. Coherent pion neutrinoproduction 



Fig. [7| shows the repeated calculations in Ref. 48] for CC 7r + production from l2 C. Only 
the diagrams with the A in s and u channels are included. We use the iV «-» A transition 
form factors in Ref. |48| to extrapolate our calculation to E v ^ 0.5 GeV 7 . The A self-energy 



modification and the pion optical potential are also the same as in Ref. 48]. This plot 
shows three different calculations. The "A unmodified, PW, 1/2 cr" calculation does not 
apply medium-modification to the A self-energy; it treats the pion wave function as a plane 
wave; and scales the total cross section by 0.5. In the "A modified, PW, 1/2 cr", medium- 
modification for the A is included. Finally, the "A modified, DW" calculation includes both 
the medium-modification and a distorted pion wave function. A good agreement between 
these results and those in Ref. [48[ is achieved, which is a justification for our numerical 
calculation. 

Now let's turn to our results for the total cross section of CC 7r + production in 

(anti)neutrino- 12 C scattering as shown in Fig. |HJ Here we make use of the meson-dominance 
form factors that are discussed in Ref. 8 . All the calculations include diagrams up to 
v — 2 |2j, with different (r s , r v ). We also show the NUANCE output for coherent pion 
production, which is scaled by 1/4. (In NUANCE, no wave function distortion is applied 
for the pion while pion absorption and rescattering are included in the subsequent step of 
NUANCE code 67J]; neutrino-induced and antineutrino-induced coherent pion production 
have the same cross section.) In both 7r + and ir~ production, the (1, 0.9) prediction is bigger 
than the (1, 1) and the (1, 1) bigger than the (0, 0). By comparing the differences among 
the three calculations with those in the incoherent productions [3j, we see the coherent 
processes are more sensitive to (r s and r v ) than the incoherent. This is consistent with the 
discussion in Sec. Ill I Al Moreover, our results are much smaller than the NUANCE output 



Ref. [48| labels form factors as C^ 2 3 and Cf 2 3 a- 
8 With the meson-dominance form factors, the conservation of the Electromagnetic currant is automatically 
satisfied in the free nucleon scattering calculation and in the coherent production calculation by using the 
optimal approximation. 
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FIG. 6: (Color online). The angle 6^ dependence of da/dCl^ of coherent ir° photoproduction from 
12 C with the photon energy fixed at E~ = 0.173, 0.235, and 0.29 GeV. The explanation for different 
calculations can be found in the text. The data are from Refs. 44l. |45|. 
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FIG. 7: (Color online). The total cross section for coherent CC 7r + production in neutrino- 12 C 
scattering. The calculation ingredients are the same as those in Ref. 48J. See the text for the 
explanation of different curves. 
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FIG. 8: (Color online). The total cross sections for coherent CC tt + (tt ) production in 
(anti)neutrino- 12 C scattering. The explanation for different curves can be found in the text. 
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FIG. 9: (Color online). The total cross section for coherent NC 7r° production in both neutrino- 
and antineutrino- 12 C scatterings. The explanation for different curves can be found in the text. 



in the two plots 9 . As we know, the previous calculations 60J, implemented by NUANCE, 
give bigger cross sections for coherent pion productions than the measured ones l47t . It is 



noticed in the Fig. [7] and Refs. 48-54, 56[ that including the medium-modification on the 
A and the distortion of the pion wave function reduces the cross section significantly. 

Fig. [9] shows our results for NC tt° production in neutrino- and antineutrino- 12 C scat- 
terings. Three different calculations are presented in the same as way as in Fig. |HJ The 
systematics in them are the same as in the CC productions: first, the (1, 0.9) gives the 
biggest cross section and the (0, 0) gives the smallest; second, cross sections are sensitive to 
r s and r v , compared to incoherent NC productions in Ref. [3J; finally our results are much 
smaller than the NUANCE output even after including the pion absorption (15% reduction). 

In addition, as mentioned in Ref. j60[, the coherent production is dominated by the for- 
ward production at the high energy region. By using the conservation of vector current (the 
leptonic current should be proportional to the momentum transfer in the forward kinemat- 
ics), the contribution of the vector current in the full hadronic current is small [see Eq. (j3J)], 



In NUANCE, an over all 15% reduction is expected after the pion absorption is included [67 1. 
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and hence the interference between the vector current and axial current is small. As the 
result, neutrino- induced and antineutrino- induced production should have similar cross sec- 
tions (see the NUANCE output shown in the plots). However this is clearly violated in 
the energy region of this paper (see results in Fig. |H] for the CC production and in Fig. M 
for the NC production). Furthermore, for an isoscalar nucleus like 12 C, the axial current 
in the hadronic CC and in the hadronic NC should have the same strength by using the 
Wigner-Eckart theorem. In the leptonic current, this ratio is \/2 : 1. Because axial current 
dominates in both CC and NC production at the high energy region, the ratio for the cross 
sections between them should be 2 : 1 in both neutrino and antineutrino scatterings (ignore 
the u and d quark mixing) (63[, which is also represented by the NUANCE output in the 
plots. But this ratio is not satisfied at low energy, if we compare the (anti)neutrino results 
in Fig. [8] with the (anti)neutrino results in Fig. [9j So, it will be interesting to extrapolate 
our low-energy results to high energy and find out the transition region where the predicted 
high-energy behavior starts to emerge. 



C. Coherent NC photon production 

In Fig. dm we show our results for coherent NC photon production in both neutrino- 
and antineutrino- 12 C scatterings 10 . In accordance with the low detection efficiency for 
low energy photons in the MiniBooNE experiment, we require the photon energy in the 
laboratory frame to be bigger than 0.15 GeV in the calculation, which on the other hand 
simplifies the calculation because of the absence of the infrared singularity. The labeling of 
curves is the same as in pion production. Here all the necessary diagrams up to v = 3 are 
included. As we know from Refs. 0, Hj], all the v = 2 contact diagrams do not contribute to 
NC photon production, and the v = 3 diagrams are due to c\ and e\ coupling mentioned in 
the previous pion photoproduction calculation (the e\ contribution vanishes for an isoscalar 
target). We can observe the effect of c\ coupling, by comparing the "(1, 1), v — 3" curve with 
the "(1 , 1), v — 3, (ci = 3)" curve [c\ = 1.5 in the calculations without {c\ = 3) labeling]. It 
increases the total cross section by roughly 10%. So, at low energy the contributions of the 
contact terms are not negligible, and a similar observation is made in photoproduction of 
pions as shown in Fig. |5j But the possibility of introducing such couplings to explain the low 
reconstructed energy excess events in MiniBooNE is not quite promising at least considering 
only low- energy neutrino contributions MJTdj. (In Ref. |5j], the number of excess events at 
low reconstructed neutrino energy is roughly two times bigger than the number of the A 
radiative decay estimated in the MiniBooNE's background analysis.) The contributions of 
these terms at high energy region still need to be studied. Moreover, the hierarchy among 
cross sections using (1, 0.9), (1, 1), and (0, 0) (with c\ = 1.5) is also consistent with the 
discussion for the Fig. However the difference among them is less significant than that 
in pion production, which is probably due to the absence of the distortion of photon wave 
function. The spreading between "(1 , 1), v — 3, (ci = 3)" and "(1 , 0.9), v = 3" gives a sense 
of uncertainty of these calculations. Furthermore, the cross section in neutrino scattering is 



10 



In NUANCE, there is no manifest coherent photon production channel. But in the MiniBooNE's analysis 



62j, the total photon production is computed by scaling the total pion production, set from the total 



NC 7T° production data, by a proper branching ratio, which in principle has the contribution from the 
coherent production. 
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FIG. 10: (Color online). The total cross section for coherent NC photon production in both 
neutrino- and antineutrino- 12 C scatterings. 0.15 Gev is applied. The explanation for 

different curves can be found in the text. 

bigger than that in antineutrino scattering, which is different from the expectation about 
them at high energy region. 

IV. CONCLUSION 

In this paper, we have studied coherent neutrinoproduction of photons and pions with 
E v ^ 0.5 GeV. This paper, combined with those in Refs. 0, S] about the productions from 
free nucleons and the incoherent productions from nuclei, complete the study on these pro- 
cesses at low energy. The series is motivated by the low reconstructed energy excess events 
in the MiniBooNE experiment. The QHD EFT (with the A introduced) has been used 
in these works. It is a Lorentz-covariant, meson-baryon EFT with a nonlinear realization 
of the chiral symmetry. The U{1)em gauge symmetry and chiral symmetry guarantee the 
conservation of vector current and partial conservation of axial current. These constraints 
seem trivial at the nucleon level, but important in many-body calculations. For example, 
various procedures would have to be applied by hand to make sure the vector current con- 
served, if gauge symmetry is not manifest. Even worse, this procedure can be entangled 



(0,0), v=3 
(1,1), v=3 




(0,0), v=3 
(1,1), v=3 
(1,0.9), v=3 
(1,1), v=3 (c 1= 3.0) 
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with the specific approximation scheme. Another advantage of working in the EFT is the 
power-counting of diagrams, through which we can address the relevance of some interac- 
tion vertices. In incoherent NC photon production, we see the two contact terms c\ and 
ex, which can be partially related with the newly proposed meson's anomalous interactions, 
are negligible (they are at next-to-next-to- leading order). Their contributions do show up 
in the coherent productions, e.g. coherent pion photoproduction and NC photon produc- 
tion, as demonstrated in this paper, but do not seem to increase the photon production as 
substantially as needed to explain the excess in the MiniBooNE experiment. 

After discussing this paper in a big context, let's proceed to summarize the specifics. The 
so-called optimal approximation is introduced to simplify the calculation of nuclear matrix 
element, in which the one-body current matrix elements are factorized out. Meanwhile the 
modification on the one-body interaction amplitude is taken into account. The real part 
of the nucleon and the A self-energies is calculated by using the mean-field approximation 
of this model. The change of the A width is parameterized in a phenomenological way 
according to pion-nucleus scattering data. The medium-modifications have been tested in 
incoherent pion production in Ref. jjjj]. The Eikonal approximation is used to handle the 
distortion of the final pion wave function. Moreover, we have compared our approximation 
with the one used in [35j in which the authors introduce other densities besides the baryon 
density used in our approximation. It is shown in Fig. [3] that the two methods give similar 
results. 

We calculate the differential cross sections for pion photoproduction, which serves as 
the benchmark for our approximations, and then calculate the total cross sections for pion 
neutrinoproductions. The results are sensitive to r s , r v , and the contact term c\. The 
disagreement at high energy with the fixed pion angle shown in Fig. and at the big 
pion angle with fixed photon energy shown in Fig. [6] seems to indicate that it is necessary 
to go beyond the one-body current approximation to explain the full data. However, to 
resolve the disagreement, both the A dynamics and the distortion of the pion wave function 
should be understood better as well. In addition, we also compare our neutrinoproduction 
results with those in the literature to check our numerical calculation. Finally in photon 
neutrinoproduction, the total cross sections also depend on r s , r v , and c\. Changing c\ from 
1.5 to 3 increases both neutrino- and antineutrino-induced photon production by roughly 
10%. 

Now, let's come back to the question about the photon production being the excess 
events in the MiniBooNE experiment. One tricky point should be pointed out here. The 
reconstructed neutrino energy is based on CC quasi elastic scattering kinematics, which 
can underestimate the neutrino energy in the photon production. So the high energy neu- 
trino contribution to the photon production should be addressed before drawing a definite 
conclusion for this question. The calculations in Ref. [3} and this paper illustrate the ap- 
proximations used in both incoherent and coherent productions, and provide important 
calibration for the modification of the one-nucleon interaction amplitude in nuclei. The 
sensible extrapolations of current results to the high energy region will be pursued in future 
work. 
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